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This paper presents a study of a minimum time-to-climb trajectory analysis for aircraft flying in a two- 
dimensional altitude dependent wind field. The time optimal control problem possesses a singular control 
structure when the lift coefficient is taken as a control variable. A singular arc analysis is performed to obtain 
an optimal control solution on the singular arc. Using a time-scale separation with the flight path angle treated 
as a fast state, the dimensionality of the optimal control solution is reduced by eliminating the lift coefficient 
control. A further singular arc analysis is used to decompose the original optimal control solution into the flight 
path angle solution and a trajectory solution as a function of the airspeed and altitude. The optimal control 
solutions for the initial and final climb segments are computed using a shooting method with known starting 
values on the singular arc. The numerical results of the shooting method show that the optimal flight path 
angle on the initial and final climb segments are constant The analytical approach provides a rapid means for 
analyzing a time optimal trajectory for aircraft performance. 


L Introduction 

The climb performance of an aircraft is an important design requirement for establishing trajectories to reach a 
specified altitude and airspeed after takeoff in some optimal manner. For transport aircraft, a climb segment may 
follow a trajectory designed to achieve an optimal fuel consumption or a minimum time. Trajectory optimization 
problems to minimize aircraft fuel consumption or time of climb had been studied by various contributors in the 7(Ts 
and SG’s. 1 ” 7 In recent years, minimum time problems have also been examined for many space systems. 8 " 10 Many 
approaches are found in literature for analyzing optimal flight trajectories for a minimum time-to-climb of an aircraft 
One such method is based on the singular perturbation method that has been investigated for the minimum time-to- 
climb problem. 3,5 The singular perturbation method performs a time-scale separation of the fast and slow states in 
flight dynamic equations so that the dimension of the problem is reduced by the order of the fast states. If the aircraft 
is modeled as a point mass with three state variables: altitude, speed, and flight path angle, then the flight path angle is 
considered as a fast state and therefore its differential equation can be treated in a quasi-steady state approximation. 14 
This allows the flight path angle be treated as a control variable for the minimum time optimal control problem. 

Another popular method is based on the energy state approximation method; 1,4T6 which facilitates a-simple means 
for computation by combining the altitude and speed variables into the energy state variable that represents the sum 
of the kinetic energy and potential energy of the aircraft during climb. Because of the order reduction in the state 
equations, the energy state approach enjoyed a popularity in minimum time optimal control problems. The total 
energy level curves thus can be viewed as curves of sub-optimal climb path. 11 Climbs or dives along the energy level 
curves theoretically is supposed to take virtually little time. In some climb maneuvers, an aircraft often has to execute 
an energy dive to trade altitude for speed. Once a desired speed is attained, the aircraft climbs out of the current energy 
level curve and transitions into a curve that takes the aircraft to the next energy level curve. This transition curve cuts 
across the energy level curves in an optimal manner such that the altitude can be gained in a minimum time. This 
curve is known as the energy climb path. Fig. 1 illustrates an energy climb maneuver. 

The energy climb path turns out to be a singular arc control problem that can be analyzed by the Pontryagin’s 
minimum principle. 12 In brief, the minimum principle introduces the concept of Hamiltonian functions in analytical 
dynamics that must be minimized (or maximized) during a trajectory optimization. The optimization problems are 
formulated using calculus of variations to determine a set of optimality conditions for a set of adjoint variables that' 
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provide the sensitivity for the control problems, 13 This is known as an indirect optimization method which typically 
involves a high degree of analytical complexity due to the introduction of the adjoint variables that effectively doubles 
the number of state variables. Furthermore, the solution method frequently involves solving a two-point boundary 
value problem. Nonetheless, the adjoint method provides a great degree of mathematical elegance that can reveal the 
structure of a problem. In some cases, exact analytical optimal solutions can be obtained. Except for simple problems, 
many trajectory optimization problems require numerical methods that can solve two-point boundary value problems 
such as gradient-based methods, shooting methods, dynamic programming, etc. 



Fig. 1 - Time- Optimal Energy Climb Path 

Within the framework of the Pontryagin’s minimum principle, the singular-arc optimal control method is an inter- 
mediate method for the trajectory optimization. The existence of a singular arc in the time optimal control can simplify 
the trajectory optimization significantly. Briefly, the singular arc is described by a switching function that minimizes 
a Hamiltonian function when the Hamiltonian function is linear with respect to a control variable. 

In this study, we will examine an aspect of the minimum time-to-climb problem for an aircraft flying in the 
presence of a two-dimensional atmospheric wind field. An analytical solution for the singular arc is obtained. Wind 
patterns at a local airport can affect the climb performance of aircraft. While the time-optimal climb problems have 
been thoroughly studied in flight mechanics, the effect of winds are usually not included in these studies. A solution 
method of a minimum-time to climb will then be presented for computing a minimum time-to-climb flight trajectory. 


IL Singular Arc Optimal Control 


In our minimum time-to-climb problem, the aircraft is modeled as a point mass and the flight trajectory is strictly 
confined in a vertical plane on a non-rotating, flat earth. The change in mass of the aircraft is neglected and the engine 
thrust vector is assumed to point in the direction of the aircraft velocity vector. In addition, the aircraft is assumed to 
fly in an atmospheric wind field comprising of both horizontal and vertical components that are altitude-dependent. 
The horizontal wind component normally comprises a longitudinal and lateral component. We assume that the aircraft 
motion is symmetric so that the lateral wind component is not included. Thus, the pertinent equations of motion for 
the problem are defined in its the state variable form as 


h = v sin 7 -f Wh 


( 1 ) 


v = 


7 


T - D — W sin 7 


m 

L — W cos 7 
mv 


4 - 


— w x cos 7 — Wh sin 7 

w x sin 7 — Wh cos 7 
u 


( 2 ) 

0 ) 
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where h is the altitude, v is the speed, 7 is the flight path angle, T i-s the thrust force, D = CpqS is the drag force, 

L = ClqS is the lift force, W is the aircraft weight, m is the aircraft mass, and w x == w x (h) and vjh = Wh ( h ) are 

the respective temporal average horizontal and vertical wind field components as functions of the altitude. Thus, the 
time rate of change of the wind field can be computed as 

w x = w'Ji == v/ x (v sin 7 4- Wh) (4) 

w h = w h h = w h (v sin 7 4- w h ) (5) 


where the prime denotes the derivative with respect to the altitude h. The gradient of the wind velocity with respect to 
the altitude is also called a vertical wind shear. 16 

The problem is now posed as a minimization of the time-to-climb from an initial speed and an initial altitude to 
a final speed and a final altitude at a final time tf. To formulate the time-optimal control problem, we consider the 
following cost function 

f tf 

J= dt (6) 


subject to state constraints by Eqs. (l)-(3). 

The boundary conditions for the problem are the initial and final altitude and airspeed of the aircraft as 


h (0) = h$ h (tf) — hf 

V (0) = UQ v (tf) = Vf 


M (h (0) , v (0)) = M 0 


M (h (t f ) , v (t f )) = M f 


where M is the Maeh number as a function of the altitude and airspeed. The flight path angle 7 (0) or -y '(tf) may be 
free or fixed. 

To solve for the time-optimal control problem of Eq. (6), we apply the Pontryagin’s minimum principle which can 
be stated as follows: 

Let a: (t) : (0, t) R n be state variables and u (t) : (0, tf) -r E m be in a set of admissible control U that 
guides a dynamical system described by x = / ( x (t) , u (t)) from an initial state x (0) to a final state x (tf) where 
/ (a:, u) : E n x E m — > R n is some function, there exist a set of variables A (t) : (0, t) E n , called adjoint variables, 
such that we have the following necessary conditions for optimality casted in a Hamiltonian canonical structure 


dH (x (t ) , u (t) , A (t)) 


(7) 


A = — (x (8) 

dx T 

where H (x (t ) , u ( t ) A (t)) is the Hamiltonian function of the dynamical system which, for a minimum time-to-climb 
optimal control problem, is defined as 


H (x (t ) , u (t) , A (t)) = 1 + X T (t) f (x (t) , u (t)) 


(9) 


The optimal control is given by the following necessary condition 

u* (t) — argminiJ (x* (t ) ; u (i) , A (t)) (10) 

u0A 

If the dynamical system is autonomous, that is / is not an explicit function of time, the Hamiltonian function is 
required to be zero throughout the optimal solution 6 


S (x ( i ) , u (t) , A (t)) = 0 


( 11 ) 


to 


15 


A singular control exists if the control variable u appears linearly iii the Hamiltonian function H corresponding 


d 2 H 

du 2 


~ 0 


( 12 ) 
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Let S (x, A) : R n x R n — > R m be a switching function defined by 


5 = 


dH 

du 


(13) 


Since the control u does not appear explicitly in S, it cannot be determined explicitly. However, by differentiating 
S repeatedly until the control u appears explicitly, then the control u is so determined. For the control u to be optimal, 
the Kelley’s condition must be satisfied 


(-if — (—) > 0 

[ J du{ dt™ J - 


(14) 


For this problem, the Hamiltonian function is defined as 

( T — JD — PF sin 7 

— — — - w x cos 7 - Wh sin 7 

m 

+ \ f L ~ w CQS 7 j u>x sin 7 - w h cos 7 A 
7 V mv v J 

where A&, A v , and A 7 are the adjoint variables. 

We define the specific excess thrust F and the wing loading factor n as functions of the altitude, Mach number M, 
and Reynolds number Re as 


F(h,v y C L ) = 


T[h y M{h,v)\- {C D> o [Refcv) j +KCj} g [p jh ) , M (fi, v)) S 


n ( b , v , Cr,) = 


PF 

C L g[p(fQ,M(M)]g 

PF 


(16) 

(17) 


where Cd,o is the profile drag coefficient and K is the induced drag parameter. 

Then, the necessary conditions for optimality result in the following adjoint differential equations 


A h 


dH 

dh 


Av 


A 0 


dF 
1 dh 

gdn 


(u) x cos 7 Fw h sin 7^ (v sin 7 ~r Wh) — (w x cos 7 -b w k sin7j w h 


v dh 


+ {w x si n7 


h cos ?) (si 


Wh\ 

sin 7 -| j + 

v / 


w ’h ( w x sill 7-w'hCOsy) 


(18) 


: _ dH _ ^ . 

A<y — _ Ah sin 7 

dv 


dF ( • ./ . \ . 

p— - cos 7 + io A sm sin 7 


A^ 


g dn gn gcps 7 (u4 sin 7^ cos 7) 


oh 


u u 2 


(19) 


A 7 = - 


= — AfcU cos 7 — A v j*—p cos 7 — v cos 27 4- w k sin 27^ + sin 7 — w h cos 7^ 

(w' x cos 7 4- sin 7^ 


A>y 


psm7 

v 


-F w x sin 27 — w h cos 27 + 


( 20 ) 


In order to determine the extremal control to achieve a fastest climb path, we take the lift coefficient Cl as a control 
variable. Then the optimal control is one that renders the Hamiltonian stationary : . .. 7- 


dH x 2 KCtqS x qS n _ A 7 

ST- = - A «- " - 0 =*> C* L = —V 

oCl m mv 2 KvX v 

4 of 16 


(21) 


American Institute of Aeronautics and Astronautics 



subject to the inequality constraint Cl, min < < Ch,max- 

The optimal lift coefficient in Eq. (21) may be solved by a two-point boundary value problem using a numerical 
method such as a shooting method or a gradient descent method- Such a numerical solution often does not reveal a 
structure of the optimal control analytically. In most cases, it is found that this problem can be approximated as a 
singular optimal control problem by malting an assumption that the induced drag parameter K is usually small and 
therefore can be neglected- Ln fact, the ideal induced drag parameter for an elliptical wing loading is given by 



( 22 ) 


where AR is the aspect ratio. For a typical transport aircraft, the wing aspect ratio is about 7 so that the ideal K 
parameter is about 0.045. Thus, this assumption is reasonable. 

Under this assumption, we see that d 2 H/dC} j — 0 and so there exists a singular control with a switching function 


S - A, 


1 mv 


Since qS/mv > 0 , we then obtain a bang-bang extremal control as 


Cl, max A 7 < 0 

Cr. — Cl, S ing A 7 = 0 

,min A 7 A” 0 


(23) 


(24) 


The bang-bang control law is thus viewed as a sub-optimal solution to the minimum time optimal control problem. 
Equation (24) states that the minimum-time-to-climb trajectory is approximated by three sub-optimal arcs. The first 
arc is a trajectory on which the aircraft flies at some initial altitude and airspeed at a maximum lift coefficient until 
it intercepts with a singular arc trajectory defined by the second switching condition. The singular arc is an optimal 
flight trajectory for the fastest climb 8 and is also called an energy climb path (ECP) since it is a path that crosses a set 
of level curves of constant energy heights E = h + v 2 /2g as illustrated in Fig. 1 . At some point on this trajectory, the 
aircraft climbs out of the singular arc and flies along a final arc with a minimum lift coefficient to arrive at the final 
altitude and airspeed. 

To find a singular control, we use the fact that 5 = 0 and S = 0 on the singular arc to establish that A 7 = 0 and 
A 7 = 0. Hence, this allows us to eliminate A 7 in Eq. (20) by solving for A/* 


A h — A v 


g ^ w' x cos 27 4 w h sin 27 Wh ( w * tan 7 w k) 
v ' cos 7 


v 


The remaining adjoint equations (18) and (19) now become 
' dF w 


Afr — A v < —g 


dh v 


+ 


/ / » " . ( *\ l /. Wh\ 

w x w h tan 7 4 w x v cos 7 + w h v s in 7 - [w h J j^sm 7 4- — J 


A-y A-y 


dF sin 7 


4- sin 


in 7 (w' x tan 7 - w h ^j ^sin 7 + 


(25) 


(26) 


(27) 


dv 1 v 

Differentiating Eq. (25), which is equivalent to computing 5 = A 7 , and then substituting Eq. (27) into the resulting 
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expression yield 


• I gF r / r /\ 1 dF g w' cos 27 4 - w' h sin 27 Wh 

A* = A V { ^ -<? + (w x tm 7 -toJt% -g— - + -* L — ^ ' V ~ ■ — 

) V A L \ V J OT V cos 7 V 


g n t , 9 \ / 

4- — — sm 7 (3 -j- tan 7 j 4- zrc/^ cos 7 - 


t?cos ^7 v\ v cos 7 


2 (^w x j tan 7 / _ Wfl x w x w h ( . w h \ f ’\ 2 f ■ , w h\ 

f — sm 7 H +■■ ■ - ■ 3sm7 + 2 — - Wh sm7H 

cos 7 V V / cos 7 V v J \ J \ v J 


COS 7 \ V J COS 7 A v J \ / j \ V / 

w" cos 27 -f wi sin 27 Wh ( w x tan J ~~ w h) / Wh \ \ /OON 

4- — sin - 7 4- • — > (28) 

cos 7 v \ v / 

By equating Eq. (28) to Eq. (26), the singular arc for a minimum tlme-to-climb is now obtained as 

dF v 2 dF 

f (h } v* 7 > n ) = (1 -ai)F+ (1 4-a 2 ) u— - — 4- a 3 n 4- a 4 === 0 (29) 

ov g on 

with the coefficients aa, a 2 , a 3 , and a 4 described by the following functions of the wind field parameters 

w h («4 tan 7 -w*) 


! (w' x cos 27 4 - n) h sin 27 ^ Wh (w' x tan 7 - w h 


g cos 7 


V * t 0 \ / W-Wh 

a 3 =•- ru r sin 7 (3 4- tarn 7 ) - 2w h cos 7 4 0 - 

g v ' . v cos 2 7 


/usin 7 4 -tn/i \ 2 jf^tan 7 wA , t „\ v ( * * 

« 4 =( j -+ v ^ tan 7-^J -- 3 Wx tan 7 - Wft + 


2 w x wh 
v cos 7 


Equation (29) is a partial differential equation in terms of the specific excess thrust F that describes an optimal 
climb path for a minimum time-to-climb solution for an aircraft flying in the presence of an altitude dependent at- 
mospheric wind field. Examining Eq. (29) reveals that there is a high degree of cross coupling between the horizontal 
wind field and the vertical wind field up to the second derivatives of the wind fields. Thus, not only, the wind field 
gradients affect the optimal climb path, but the wind field curvatures also play a role as well. Equation (29) results 
in a parabolic equation in terms of Cl , which can be solved to give a feedback control for the lift coefficient on the 
singular arc as a function of the three state variables h , v , and 7 . 

A certain simplification of Eq. (29) can be made by considering the concept of fast and slow states in flight 
dynamics. Ardema had shown that in optimal trajectory analysis, the three-state model of a point mass aircraft exhibits 
a time-scale separation behavior whereby the state variables h and v possess slower dynamics than the state variable 
7. 3 This time scale separation is normally treated by a singular perturbation analysis by replacing the fast state equation 
with the following equation 

L — PE cos 7 w x sin 7 — uth cos 7 


where e is a small parameter. 

Equation (34) has an inner solution and an outer solution, which is the steady state solution obtained by setting 
e = 0. The inner solution can usually be solved by the method of matched asymptotic expansion. The inneir solution 
is also called a boundary -layer solution in reference to the historical origin of the singular perturbation method in 
the fluid boundary layer theory. The overall solution is dominated by the outer solution with the inner solution only 
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affecting a small initial time period. As a first order approximation, the inner solution can sometimes be ignored. In 
this case, Eq. ( 34 ) results in the following load factor 

{v sin 7 4 Wh) ( w x sin 7 — w h cos 7 ] 

n — cos 7 — — • ■■■ ■ — — ( 35 ) 

9 

Equation ( 35 ) thus dethrones the number of state equations from three to two by converting the state variable 7 
into a control variable in place of the lift coefficient Cl- Substituting Eq. ( 35 ) into Eq. ( 29 ) now yields the optimal 
climb path function 


f(h,v,j) = F 


w h ( w 'x tan 7 — ^4 ] ] dF I v, cos 27 + sin 27 ) w h (w x tan7 - w h ) 
1 — — +u— - 1 + — — - 


v^dF _ v 
9 dh g 


dv 


r 9 / VJ T Wh 

2ul. sin 7 tan 7 4- w h cos 2 j 4 — 

x n v cos 7 


g cos 7 


4 - v 


v sin 7 4- Wh 


(«C tan 7 - 


^ v ( v sin 7 4 - Wh 


(w' x tan 7 — Wfrj | — tan 7 cos 27 4 - 2 w h cos 2 7 4 


v cos 7 


— 0 ( 36 ) 


We have thus reduced the optimal climb path by removing the dependency on the lift coefficient. We now examine 
some special cases of the general optimal climb path function f (k, v, 7): 


1 , 


Steady Wind Field: 

For a steady wind field, the gradients and curvatures vanish, thereby reducing the optimal climb path function 
to the following 


/ (h, v) = F -h v 


dF 

dv 


9 dh 


0 


( 37 ) 


Equation ( 37 ) is the well-known result for the ECP without an atmospheric wind field. Thus, the optimal climb 
path in the presence of a steady wind field is effectively the same as that without the wind field effect. This result 
is not surprising and can be explained by the fact that since the inertial reference frame is attached to the air 
mass, whether the air mass is moving at a constant speed or remains stationary, the speed of the aircraft relative 
to the air mass is the same, thus resulting in the same optimal climb path. 

2. Horizontal Wind Field: 

In the presence of a horizontal wind field only, the optimal climb path function becomes 


/ (7i, v, 7 ) = F 4 v 


dF 

dv 


14 


vw x COS 27 
g cos 7 


v 2 dF 2 vw' x sin 2 7 tan 7 v 3 w x sin 2 7 tan 7 
9 dh g + ~ g 2 


V 2 (4) 


sin 7 tan 2 7 cos 27 == 0 ( 38 ) 


Specializing Eq. ( 38 ) for small flight path angles by invoking sin 7 7 and cos 7 ** 1 yields 


f(Kv) =F + v 


d£ 

dv 



iP&F 
9 dh ‘ 


2 vw' x 

9 



= 0 


( 39 ) 


3 . Vertical Wind Field: 
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In the presence of a vertical wind field only, Eq. (36) becomes 


/(M, 7 ) = F 1 + 




dv 


^ v dF ( ^ ^ vw h sin 2j t WhW h ^ v 2 dF vw h cos 27 


g cos 7 


5 / 9 dh 


// ( v sin 7 4- w/t 


v 2 2v \ 


2 

1 cos 2 7 

1 

( usin 7 + Wh\ 

) 9 

l 9 ) 


= 0 (40) 


For small flight path angles, Eq. (40) reduces to 


f(h, 


i + «i N | + »K + + 


T r * 


9 ) g dh 


2v 


(»;)' 


vj 4- Wh 


0 (41) 


Vertical wind field is especially important for microburst problems. A microburst is a wind shear disturbance 
characterized by an unusually strong downdraft near the ground surface that presents an extreme hazard to 
aircraft during take-off or landing. 


in. Optimal Solution 

During a climb, the aircraft flies under a maximum continuous thrust from take-off to a point in the flight trajectory 
envelope where it intersects with the optimal climb path. The aircraft then maintains its course along the optimal climb 
path until it reaches some point on the optimal climb path where it climbs out to the final altitude and airspeed. Thus, 
there exist three climb segments during a climb as illustrated in Fig. 1 labeled as A, B, and C where B is the singular 
arc optimal climb path. The suboptimal solutions for the climb segments A and C may be defined by lines of constant 
energy heights E = b + v 2 /2 g. n 

From the optimal control perspective, the initial and final climb segments are to be determined by requiring the 
Hamiltonian function as defined in Eq. (15) to be zero. In addition, the adjoint variable A 7 is no longer restricted 
to zero according to the first and last switching conditions in Eq. (24). Thus, in general, the optimal solutions for 
the initial and final flight segments are considerably more complex than the optimal climb path solution and usually 
involve solving a two-point boundary value problem. 

First, we shall consider the singular arc optimal solution. Along this singular arc, all the state and adjoint variables 
are functions of the altitude h, airspeed v, and flight path angle 7 . However, it can easily be shown that the flight path 
angle 7 on the singular arc optimal energy climb path is in turn a function of the altitude h and airspeed v. Therefore, 
the two variables h and v uniquely determine the optimal climb path. In particular, for the case of small flight path 
angles for which Eqs. (39) and (41) apply. The optimal climb function / (h, u, 7 ) can be written as 

/(M. 7) 

where fi are explicit functions of v and A. 

The time derivative of f (h,v i 7 ) is evaluated as 

3 

/=E 

* t =0 

taking into account that the assumption 7 — 0 is built into the climb path function. 

Substituting Eqs. (l)-(3) into Eq. (43) then yields the following polynomial equation that can be solved for, the 
flight path angle 7 ' 

. 3 

C 57 5 + C 47 4 -f a (h r v) 7 * == 0 (44) 

i— 0 
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where the coefficients c* are 


( 45 ) 


(h \ ' 

<*(M) = —fo vw h 


S 6 9/3 / ■/ / \ 9/2 / 

c 4 ( A > 5 ~ “sr* 


dh 

dh 


Ci ( h , v ) 


5/i , dfi 

dh Whnr dv 


f( p »- 


ft, («k +»»»! + »)- ^ 

, , 0/*-i 9/i-.i f / , ' , ) 

WhW x ] + ~mr v “ “^7“ \ vt£? * + + pj 


dfi—2 


~vw h 


(46) 


(47) 


9/i 9v V x " A 9t; 

Equation (44) gives the flight path angle as a function of the altitude and speed. Thus the optimal climb surface 
function / (h t v y 7 ) can be replaced by a climb path function / (h y v y 7 (h, v)) upon embedding the solution of the 
flight path angle from Eq. (44): 

The adjoint variables along the optimal climb path are now determined from enforcing the conditions H ~ 0 and 


A 7 .= 0 on the singular arc as 


X h = 


9 4 vw x 4 2vw' h y - w h w x y + Whto h 

gv ( F - 7 ) 4 (v 7 4 w h ) (. g 4- vw' h 7 - w h w' x 7 4 


— 


(48) 


(49) 


gv (F - 7 ) + (vy 4 Wh) (g 4 ^7 — Whto^y 4 

The foregoing analysis has shown that the along the singular arc, the time optimal trajectory is known in the v — h 
plane as well as the flight path angle and the adjoint variables. This information greatly facilitates the optimal solutions 
for the climb segments on the initial and final arcs. Since the initial conditions on the initial climb segment are known 
according to the problem statement and the fact that its solution must terminate on the optimal climb path function 
for which all the state and adjoint variables are known, then the optimal solution for the initial arc can be computed 
using a shooting technique to integrate backward from some point B on the optimal climb path to the initial point 
A as shown in Fig. 1 , Likewise, to compute the optimal solution for the final arc, a shooting method may be used 
to integrate forward from some point C on the optimal climb path to the final point D that terminates at the desired 
altitude and airspeed. A shooting method may be established as follows: 

The flight path angle 7 not on the optimal climb path function can be solved by setting the Hamiltonian in Eq. (15) 
to zero with the usual small angle assumption, thus resulting in the following quadratic equation 


X v vw h j 2 + \ h v - A„ (g + vw x - w h w'hj 7 + 1 + ^hWh + K (Fg - WhW^ = 0 


(50) 


Equation (50) is then used to eliminate the flight path angle expression in the adjoint equation in the adjoint 
differential equations (18) and (19), which can then be parameterized in terms of the altitude as the independent 
variable instead of time using the following transformation with the small angle assumption 


d\ h ~ A hw'h. - A « [?f - («£ + 4r) (n + ^h) - («4 + «yy) *4 

dh h 1*7 + Wh 

d\ 1 _ X L _ -Afe 7 - At, .[g|g - (vj + nyr) 7 ] 

dv h vy 4- Wh 

In addition, the airspeed can also be parameterized as a function the altitude as 


dv 

dh 


v 

h 


(F ~ 7) 9 
vy 4 Wh 


- w x + W h j 


(51) 


(52) 


(53) 


Equations (50)-(53) are then integrated either forward or backward using the known state and adjoint variables on 
the optimal climb path as the starting values at points B and C. The shooting method then iterates on these known 
values on the optimal climb path until the conditions at points A and D are met 
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TV. Numerical Example 


We compute the time optimal climb trajectory of a transport aircraft in a horizontal wind field. The aircraft has a 
maximum thrust T max which varies as a function of the altitude as 


T — 'Tmax 



(54) 


where p is the density, p 0 is the density at the sea level, and a is a constant 

Using Eq. (35) with the small flight path angle approximation, the specific excess thrust is computed as 

r TrnazfpY CpjpfS 

W \poJ W 

The optimal climb path function, Eq. (39), is equal to 


(55) 


f(v,h) 


T max f PY (i _ ay2 d P\ 

W [poj { pg dhj 


C D \pv 2 S 


W 


v 2 dp 2vw 0 


pg dh 


2vw x 

9 


v z wl 


+ 


v 2 (w'P) 


y 3 ^0 (56) 


The aircraft has trie following values: T max =40000 lb, W =200000 lb, S =1591 ft 2 , Ct) ~ 0.0343, and a = 0.7. 
We will compute the singular arc for two typical horizontal wind field problems. The first wind field problem is a 
low-altitude wind shear disturbance described by the following model 

_ ^ _ fe 

w x = a i.e - a2e ^ (57) 


where, for the problem, a* = 250 ft/sec, a2 — 200 ft/sec, hi = 6000 ft, and h<i = 1000 ft. 

The second wind field problem is a high-altitude wind field problem described by the following power law model 


w x — a 3 


fh~h 3 

\ '/* 


a 


(58) 


where, for the problem, a 3 = 250 ft/sec, h 3 = —1000 ft/sec, /14 = 30000 ft, and a — 
The two types of wind field are plotted in Fig. 2. 


xIO 4 



Fig. 2 - Wind Shear and Power Law Models 
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As can seen, the wind shear profile is characterized by a strong wind field gradient near the ground which can be 
hazardous to aircraft during take-off and landing. 

To find the singular arc, we decompose the optimal climb path function, Eq. (56), in the individual functions 
fi(h>v) according to Eq. (42). For a horizontal wind field problem, Eq. (44) is then formulated as a degree 
polynomial function in terms of the flight path angle 7 . A physically correct root is obtained that yields the flight 
path angle on the singular are. The singular arc function / (A, v) then becomes a nonlinear function in terms of h 
and v. To compute the singular arc trajectory in the v — h plane, we find the zero solution of this function using a 
Newton-Rap hson method as follows: 

-1 

fiKvi) (59) 

The computed singular arc time-optimal climb paths for both wind field profiles are plotted in Fig. 3. For reference, 
we also compute the singular arc for zero wind disturbance. As can be seen, the singular arc climb paths are steep 
paths that rapidly increase the altitude with a relatively smaller change in the air speed for the case of no wind. The 
power law profile follows a similar pattern as the no-wind case, although the ground speed intercept is less. This 
would mean that with the power law profile, the aircraft has to enter the singular arc climb path at lower speed than 
if there were no wind. The singular arc for the wind shear profile is the most interesting in that the speed variation is 
significant over a wide range at a very low altitude. This is due to the strong wind field gradient over a short altitude. 
At high altitude, the three singular arcs are converging, so the effect of wind field is less pronounced at high altitude. 


Vi+l = Vi - 


df (hj, Vi) 
dv 



v, ft/sec 


Fig. 3 - Singular Arc Time- Optimal Climb Paths 

Once the singular arc has been determined, the values of h and v along the singular arc are plugged into the 4? h 
degree polynomial to solve. for the flight path angle along the singular arc. The lift coefficients are also computed. 
Figs. 4 and 5 are the plots of the flight path angle and lift coefficient on the singular arc climb paths. The flight path 
angle for the no- wind case generally decreases with altitude. At high altitude, the flight path angle for the two wind 
profiles converge to that of the no-wind case. The wind shear case as usual shows a drastic change in the flight path 
angle along the singular arc at low altitude. The lift coefficient generally increases with altitude for the no- wind case. 
The lift coefficient for the wind shear case varies greatly at very low altitude and is quite large at ground level due to 
; the low ground speed intercept required to enter the climb path. 

We next compute a complete trajectory from take-off to some final altitude and air speed. We only consider the 
wind shear case. The initial ground speed at take-off is about Mach 0.2 or 224 ft/sec and the desired air speed at 
15000 ft is Mach 0.5 and 0.6. The solutions not on singular arc require solving a two-point boundary value problem. 
However, since the adjoint variables are completely determined on the singular arc, we can solve the trajectory off the 
singular arc quite easily using a shooting method by integrating Eqs. (51) to (53) either forward or backward starting 
from the singular arc. The adjoint solution on the singular arc is plotted in Fig. 6 . 
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Fig. 4 - Flight Path Afrgie on Singular Arc 



Fig. 5 - Lift Coefficient on Singular Arc 




Fig. 6 - Adjoint Solution on Singular Arc 
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The complete trajectory is plotted in Fig. 7. To climb to 15000 ft and Mach 0.5, the climb trajectory is comprised 
of three segments. The first segment, segment AB, is the take-off segment on which the aircraft flies from take-off 
ground speed to intercept the singular arc at point B. This occurs at a very low altitude. The second segment, or 
singular arc segment BC, is the minimum-time to climb path that takes the aircraft to a higher altitude in a fastest time. 
At some point on this singular arc, the aircraft begins to depart at point C and flies on the final segment, segment CD, 
to the final altitude and airspeed. We note that the three segments join together and are tangent at points B and C. On 
this path, the departure slope is to the left of the singular arc. In the region to the left of the singular arc, because, the 
altitude is continuously increasing, the thrust is maintained at the maximum value. 

The situation for Mach 0.6 corresponding to the final segment EF is different. Since the aircraft must arrive at 
Mach 0.6 which is to the right of the singular arc. Because the departure slope must be tangent to the singular arc and 
curve to the left, there is no point below the final altitude where this tangency exists. As a result, the aircraft must fly 
past the final altitude and then reduce the engine thrust to let the potential energy be convert to kinetic energy once the 
aircraft begins to slow down. For the problem, we use a first-order model to describe the engine thrust reduction from 
the maximum value to an idle value at roughly 20% of the maximum thrust. The segment EF is sub-optimal in that the 
engine thrust varies so that it becomes a control variable in addition to the flight path angle. Nonetheless, it is a very 
good approximation since the segment EF approximately follows a total energy curve. 



Fig. 7 - Climb Trajectory 
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Fig. 8 - Lift Coefficient on Climb Trajectory 
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The lift coefficient variation on the climb trajectory is shown in Fig. 8. The lift coefficient varies greatly from 
a maximum take -off value to about 0.4 at a very low altitude. Depending on the wind shear profile, this may not 
be feasible due to the slow retraction of flaps and slats deployed during take-off. The flight path angle on the climb 
trajectory is plotted in Fig. 9. The flight path angle for the initial segment AB and the final segment CD are computed 
using Eq. (50) as a function of the adjoint solution computed from the shooting method. It turns out that the computed 
flight path angle is very nearly constant as shown in Fig: 10. Thus, an approximate trajectory can easily be computed 
by only integrating Eq. (53) using a constant flight path angle at a point of tangency to the singular arc. The final 
segment EF is computed using this approach. 

To examine the optimality of the computed solution, we compute the values of the Hamiltonian function along 
various climb segments. The segments AB, BC, and CD are optimal so the values of the Hamiltonian function are 
equal to zero. The climb segment EF is sub-optimal and that fact is clearly illustrated by the non-zero value of the 
Hamiltonian function. Nonetheless, it is a reasonable approximate solution, given that the problem would have been 
formulated with the engine thrust as an additional control variable, which would result in a more complex problem. 



Fig. 9 - Flight Path Angle on Climb Trajectory 
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Fig. 10 - Flight Path Angle Computation 
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Fig. 1 1 - Values of Hamiltonian Function 


V. Conclusions 

A solution to a minimum time-to climb problem in an altitude dependent two-dimensional wind field has been 
presented. This problem possesses a singular control structure when the lift coefficient is taken as a control variable 
with the induced drag effect neglected. The optimal climb path on the singular arc is obtained as a function of the three 
flight state variables and the lift coefficient as a control. A time-scale separation is used to reduce the dimensionality 
of the optimal climb path function by converting the lift coefficient control into the flight path angle control. Using 
the singular arc analysis, this three-dimensional function can be further decomposed into a solution for the optimal 
flight path angle and a two-dimensional climb trajectory as a function of the airspeed and altitude. The known adjoint 
variables on this optimal trajectory significantly simplifies the general optimal control solution for the initial and final 
climb segments. A shooting method is formulated to solve numerically the trajectories of the initial and final segments 
using the known solutions on the optimal climb trajectory. A numerical example for a wind shear profile is computed 
to demonstrate the analytical method. It is found that the computed flight path angle off the singular arc is nearly a 
constant and therefore can be used as a very good approximation. The analytical approach provides a rapid means for 
estimating climb trajectory of aircraft flying in the presence of atmospheric wind. 
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